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We examine the matrix theory representation of DO-brane dynamics 
at finite temperature. In this case, violation of supersymmetry by 
temperature leads to a non-trivial static potential between DO-branes 
at any finite temperature. We compute the static potential in the 1- 
loop approximation and show that it is short-ranged and attractive. 
We compare the result with the computations in superstring theory. 
We show that thermal states of DO-branes can be reproduced by matrix 
theory only when certain care is taken in integration over the moduli 
space of classical solutions in compactified time. 
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Figure 1: Anmilus diagram for DO-brane interactions. The bold lines represent world-lines of the DO-branes, 
separated by the distance L, which go along the periodic temporal direction. They bound the string world- 
sheet. 



1 Introduction 

Dirichlet p-branes are p+1- dimensional hypersurfaces on which superstrings can begin and 
end (see0i for a review). The low energy dynamics of an ensemble of N parallel Dp-branes can 
be described by the U(N) supersymmetric gauge theory obtained by dimensional reduction 
of ten dimensional supersymmetric Yang-Mills theory to the p+1- dimensional world- volume 
of the branei The Yang-Mills theory gives an accurate perturbative representation of the 
Dp-brane dynamics when the separations between the branes is largeBBi It represents a 
truncation of the full string spectrum to the lowest energy modes. The full string theoretical 
interactions between a pair of Dp-branes is computed by considering the annulus diagram 
shown in fig. [I]. The short distance asymptotics of this diagram are dominated by the open 
string sector whose lowest modes are the fields of ten dimensional supersymmetric Yang-Mills 
theory. On the other hand, long distance asymptotics are most conveniently described by the 
dual description of this diagram as a closed string exchange and the relevant field theoretical 
modes are those of ten dimensional supergravity. That these are also represented by the 
dimensionally reduced super Yang-Mills theory is a result of supersymmetry and the fact 
that, for fixed Dp-brane positions, the ground state is a BPS state. At zero temperature, 
because of supersymmetry, the interaction potential between a pair of static DO-branes 
vanishes independently of their separation. Their effective action has been computed in an 
expansion in their velocities, divided by powers of the separation and is known to bell 

This result agrees with the effective potential for the interaction of DO-branes in ten dimen- 
sional supergravity. Note that, for weak string coupling, the DO-brane is very heavy. 

In this paper, we shall consider the description of DO-brane interactions in type IIA super- 
string theory using matrices. Even at very low temperatures, non-BPS states are important 
to the leading temperature dependence. We perform 1-loop computation of the effective in- 
teraction between static DO-branes in the matrix theory at finite temperature and compare 
with the known superstring computations. We show that the results of the two computa- 
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tions are similar in the low temperature limit but an extra integration over the temporal 
component of the gauge field, is present in the matrix theory. At finite temperature, because 
the Euclidean time is compact, the temporal gauge field can not be removed by a gauge 
transformation. This integration is needed in order to describe correctly thermodynamics of 
DO-branes both in the matrix and superstring theories. 

The paper is organized as follows. In section 2 we discuss the formulation of the matrix 
theory at finite temperature. In section 3 we perform one loop computation of the effective 
interaction between static DO-branes at finite temperature and show that it is attractive, 
and short-ranged. In section 4 we compare this result with the superstring computations 
and discuss the conditions under which the two computations agree. Section 5 is devoted to 
the discussion of our results and, in particular, the origin of the divergence of the classical 
thermal partition function of DO-branes which is cured by quantum statistics. 

2 Matrix theory at finite temperature 

We shall consider the matrix theory description i of the effective dynamics of DO-branes in a 
type-IIA superstring theory which is derived by the reduction of ten dimensional supersym- 
mertic Yang-Mills theory which has the action 

S YM [A, 8} = ^-J drTR Qi£ + (2) 

to zero spatial dimension: = A^ij) , 6 = 6{r). 

The thermal partition function of this theory is given by 

Z YM = / [dA(r)][de(r)} exp (S YM [A, 6}) (3) 

where Sym is the Euclidean action and the time coordinate is periodic. The bosonic and 
fermionic coordinates have periodic and anti-periodic boundary conditions, 

A li (r + P)=A ll (T), (4) 
9(t + P) = -9(t), (5) 
P = l/k B T, (6) 

where T is the temperature and hs is Boltzmann's constant. Gauge fixing will be necessary 
and will involve introducing ghost fields which will have periodic boundary conditions. 

The representation ([|) of the thermal partition function can be derived in the standard 
way starting from the known Hamiltonian of the matrix theory! and representing the thermal 
partition function 

Z YU = tie- f3H (7) 
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via the path integral. The trace is calculated here over all states obeying Gauss's law 
which is taken care by the integration over A in (|3|). This representation of the matrix 
theory at finite temperature have been already discussed000, but the temperature induced 
interaction between .DO-branes described below was never identified. 

In matrix theory, the diagonal components of the gauge fields, a a = A aa , are interpreted 
as the position coordinates of the a-th DO-brane and they should be treated as collective 
variables. Static configurations play a special role since they satisfy classical equations of 
motion with the periodic boundary conditions and dominate the path integral as g\ M —>■ 0. 
Notice that there are no such static zero modes for fermionic components since they would 
not satisfy the antiperiodic boundary conditions. [] This is an important difference from 
the zero temperature case and a manifestation of the fact that supersymmetry is explicitly 
broken at non-zero temperature. 

In the following, we will construct an effective action for these coordinates by integrating 
the off-diagonal gauge fields, the fermionic variables and the ghosts, 

S cS [a a ] = - In / [da%] J] [dAf] [d9] [rfghost] exp (S YM - S gi - S gh ) . (8) 

Generally, this integration can only be done in the a simultaneous loop expansion and expan- 
sion in the number of derivatives of the coordinates a a . Such an expansion is accurate in the 



limit where 



a? 



are large for each pair of DO-branes and where the velocities are small. 
Since these variables are periodic in Euclidean time, small velocities are only possible at low 
temperatures. The remaining dynamical problem then defines the statistical mechanics of a 
gas of DO-branes, 

2ym = / Yiid^ir)] exp (-S eG [<F}) . (9) 

We expect that the zero temperature limit of S e s reduces to ([]]). 

We shall find several subtleties with this formulation. If the effective DO-brane action is 
to reproduce the results of a string theoretical computation, the integration over must be 
performed in both cases. 

The effective action is a symmetric functional of the position variables a a (r). Only the 
configuration of these coordinates needs to be periodic. Therefore the individual position 
should be periodic up to a permutation. The variables in the path integral ® should 
therefore be periodic up to a permutation and the integral should be summed over the 
permutations. 



b This is a difference between our computation at finite temperature and computations of the Witten index 
for the matrix theory where fermions obey periodic boundary conditions. 
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3 One loop computation 

We will compute the effective action S e g in a simultaneous expansion in the number of loops 
and in powers of time derivatives of the DO-brane positions. 

We decompose the gauge field into diagonal and off-diagonal parts, 

Af = a a ^ + g YM Af (10) 



where A^, a = so that the curvature is 



5^f^ + g Y MDfAf-g YM DfAf- ig 2 YM ' T 



■pa/3 



Aft, Ay 



ir. 



where 



and 



ft, = - ®K ( 12 ) 



Df = d li -i{fi- ag) • (13) 

In the Yang-Mills term in the action, we keep all orders of the diagonal parts of the gauge 
field and expand up to second order in the off-diagonal components, 

fc) =E~r- te)+2E^ Q UX^f - D^Df + K (/" - fr)) Af+. . . 

9YM a 9YM a/3 \ ) 

(14) 

We will fix the gauge 

DfAf = 0. (15) 
This entails adding the Fadeev-Popov ghost term to the action 

S* = f E {** (-Dff </* + igvu^Df \A„ c] } (16) 

J a/3 1 ' 

There is a residual gauge invariance under the abelian transformation, 

< oJ + d,x a - (17) 
We shall use this gauge freedom to set the additional condition 

d a% = (18) 

and to fix the constant^ 

- ir//3 < a* < tt//3. (19) 

c These constants are related to the eigenvalues of the holonomy 



Pexp (1tAo(t) \ =fi t diag (e i/3a o, 



known as the Polyakov loop winding along the compact Euclidean time. It can not be made trivial by the 
gauge transformation if T ^ 0. 



The ghost for this gauge fixing condition decouples. 

Keeping terms up to quadratic order in A, c, c, 9, the action is 

S = J ( A~ E (/£,) 2 + \ E K a (-^D^Df + D**D? + 2* (/" - /£,)) Af 



+ ]T <? a {Df) 2 c aP + -9 Pa lfl D^9 a P ) . (20) 

a/3 ^ 



The effective action obtained by integrating over A, c, c, 9 is 

S cS = /E ^(/^) 2 + E {^TRln {-S^{pff + 2i(%, - $ 



— TRln {-{Dff) - -TRln (i^D?)) . (21) 



5.i Leading order in time derivatives 

We will evaluate the determinants on the right-hand-side of ( pi]) in an expansion in powers 
of the derivatives of a(r). The leading order term can be found by setting a = const.. In 
this case, /" = and (here we retain the tree-level term with time derivatives) 

S CS = 8J2 {TR B ln (-(Dff) - TR F ln (-(Dff)} (22) 

Q</3 

where the subscript B denotes contributions from the gauge fields and ghosts, whereas 
F denotes those from the adjoint fermions. The determinants should be evaluated with 
periodic boundary conditions for bosons and anti-periodic boundary conditions for fermions. 
(Note that, because of supersymmetry, if both bosons and fermions had identical boundary 
conditions the determinants would cancel. This would give the well-known result that the 
lowest energy state is a BPS state whose energy does not depend on the relative separation of 
the DO-branes.) The boundary conditions are taken into account by introducing Matsubara 
frequencies, so that 
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Using the formula 



we obtain the result 



"T + r-T (24) 



e C11 = e 



da% yj ( cosh/?|a Q - aP\ + cos/? ( 



/ ^0 TT ■ LV^ »J (25) 

' 2n a<p \coshp\a a - 3P\ - cosfi (ag -aS) 
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In order to find the effective action for a a , it is now necessary to integrate the temporal 
gauge fields over the domain (— 7r//3,7r//3]. This integration implements the projection 
onto the gauge invariant eigenstates of the matrix theory Hamiltonian. 

In the case where there is a single pair of DO-branes, N=2, the integration over in 
(p5[) can be done explicitly to obtain the effective action 



S, 



where 



P(z) = 1 + 2Alz z + 12649z 4 + 254009;?° + 24349012 8 + 12456773^ 10 + 361191812 12 
+611785892 14 + 6191459z 16 + 36109171z 18 + 124627792 20 + 2432171z 22 

+254919^ 24 + 12439^ 26 + 271z 28 - z 30 , (27) 

z = exp(— ,51a 1 — a 2 |) and we have included the tree level term, which gives the non-relativistic 
kinetic energies of the DO-branes. 

The effective action has the low temperature expansion 

S cS = f dr ( -i- V (k a Y - \ (l^e-^-^\ _ 16384e- 4 ^"l 

+ 561 f 36 e-^ 1 -^l + ...)) (28) 

We shall compare in the following section this result with the superstring computation of 
the effective interaction between DO-branes. 

4 String theoretical interactions 

The effective interactions of DO-branes in superstring theory is given by computing the 
annulus diagram shown in fig. III. This was done in ref . (and in0 for Dp-branes). The 
result of summing over physical (GSO projected) superstring states gives the free energy 



/■KOL JO 

where 



v 



02 ^ n ft±£V m 



irot'l ) £± x \ 1 - e 



G 2 (v\iz)= Yl e~^ 2q+1) (30) 

Q= — OO 

L is the brane separation and v is a parameter which weights the winding numbers of strings 
around the periodic time direction. An extra factor of 2 accounting for the exchange of the 
two ends of the superstringlll ending on each of the two DO-branes is inserted in 



The product in the integrand represents the sum over string states, with requisite degen- 
eracies, 

oo / I i -nl \ 8 oo 
n=l \ c / N=0 

where is the degeneracy of the either superstring state at level N. For the lowest few 
levels, do = 8 and Eq = L/2ira' . 

Inserting fl3~l|) in (|29|) and integrating over Z, the free energy has the form 



2 00 

F(/3,L, I /) = -X;d J v In 

^ N=0 

where the string energies are given by the formula 



\ _ e -f3E N +i-wu 



~\_ _|_ e —(3E N +i-KU 



(32) 



2na' V L 2 

This results in the partition function 



L / 4vrViV . , 



Z str (/5,L,z/) = e- /3F = 1] 

N=0 



(34) 



The physical meaning of the last formula is obvious: the partition function equals the ratio 
of the Fermi and Bose distributions with the power being twice the degeneracy of states and 
iv playing the role of a chemical potential. The factor of 2 in the exponent 2d^ in (R3j) is 
due to the interchange of the superstring ends as is already mentioned. It will provide the 
agreement with the matrix theory computation. 

In order to compare with the Yang-Mills computation, we should first re-scale the coor- 
dinates so that the mass of the DO-brane appears in the kinetic term as in (p. The mass is 
given by the formula 

M = —^= (35) 
g s Va' 

and the Yang-Mills coupling gywi is related to the string coupling g s by the equation 

9YM = 4^)3/2- ( 36 ) 

The physical coordinate of the a-th DO-brane is identified with 

g° = 2Wa Q . (37) 



Taking N=2 in (p5|) and identifying L = 2ira'\d 1 — a 2 \, we see that the integrand in (p5| 
coincides with ( |34"D truncated to the massless modes (N = 0) provided u — (3{a\ — a 2 ,). 
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It is clear that the integral over is responsible for the "mismatch" of the effective 
actions between the string theory computation and matrix theory computation of the free 
energy. In the string theory done in the spirit of ref . 0, the parameter v appears in the same 
place as /3(aJ — afy/n but is not integrated. It is associated with the interaction of the ends 
of the open string with an Abelian gauge field background A^(t,x): 

SWt = / dafAr (38) 
If the ends are separated by the distance L, e.g. along the first spatial axis, then 

u= [ dr (A>(r, 0, . . .) - A (r, L, . . .)) (39) 



o 



since x^(r) = (1,0) on the boundaries. The matrix theory automatically takes into account 
the integration over the background field while in the string theory calculation of ref. the 
background field is fixed. This is just a reflection of the fact that matrix theory is an effective 
low-energy theory of D-branes, while the older string theory did not treat the boundaries as 
dynamical objects. However, it is interesting to notice how close some of the earlier string 
papers came to such a description simply by the requirement of consistency 0. Further, in 
the context of matrix theory it is natural to take the exponential of the free energy (|32|) as 
in eq. fl34]) , and only integrete over v afterwards, a procedure not entirely obvious in a string 
theory where the boundaries are not dynamical objects. This will be discussed further in 
the next section. 

An exact coincidence between the matrix theory and superstring results is possible only 
when the higher stringy modes are suppressed. Usually, the truncation of the string spectrum 
to get Yang-Mills theory is valid for small a', that is when we are interested in temperatures 
which are much smaller than l/y/a'. In fact, the condition in our case is a little different than 
this once the length L appears as a parameter in the spectrum (0). Then, the spectrum 
can be truncated at the first level only when 



which is the energy gap between the first two levels. If the temperature is small, this condition 
is always satisfied unless the length L is not too large. In other words the truncation of the 
spectrum to the lightest modes is valid for (3 ^> L (or TL Cl). 

It is also interesting to discuss what happens in the opposite limit L ^> (3 where the 
interaction between DO-branes is mediated by the lightest closed string modes. The super- 
string free energy can be evaluated in this limit by the standard modular transformation 
which relates the annulus diagram for an open string with a cylinder diagram for a closed 
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string. Introducing the new integration variable s = 2ir 2 /l, we rewrite (|29|) as@ 



nrr „ , 87T 4 f°° ds -L 2 /2sa'^ ( 



i(3 2 s 



-(2n+l)s N 



2tt 3 



ol 



n — ■ < 41 > 



n=l 



In the limit where the brane separation is large the integration over s is concentrated in the 
region of large s ~ L 2 . Substituting the large- z asymptotics 

6 2 {y \iz ) -> 2 cos (vrz/)e-^ /4 (42) 

and evaluating the saddle-point integral, we get 

F[L, /3, i/] oc cos "f/^ e -iVw/W - (43) 

Exponentiating and integrating over z/, we have 

1 ^Z str (/3, L, i/) oc ^(^ 2 - 8vr 2 a 3 e - Lv ^Zi^7 /W- (44) 
i L 8 



Taking into account (|37| ) the exponent at the low temperatures is the same as in ( p5|) but 
the pre-exponential differs. The dependence of the pre-exponential on L in the superstring 
case emerges because the splitting between energy states in (|33D is of order 1/L and the 
truncation condition (EO) is no longer satisfied when L is large. Higher stringy modes are 
then not separated by a gap and the continuum spectrum results in the L-dependence of the 
preexponential. As usual, the limits of L — ► oo and T — > are not interchangeable in the 
superstring theory. 



5 Discussion 

Our main results concern DO-brane dynamics at finite temperatures. We have computed 
the 1-loop effective action for the interaction of static DO-branes in the matrix theory at 
finite temperature and compared it with the analogous superstring computation. We have 
seen that an extra integration over the eigenvalues of the holonomy along the compactified 
Euclidean time is present in the matrix theory. The two computations agrees in the low 
temperature limit provided the superstring thermal partition function is integrated over the 
Abelian gauge fields ao's living on DO-branes. 

The integration over a 's is of course natural in the context of the Yang-Mills theory, 
where it expresses that only gauge-invariant states should contribute to TKe~^ H . But it 
is also natural from the point of view of the DO-brane physics. It can be seen as follows. 
Suppose we make a T-duality transformation, which interchanges the Neumann and Dirichlet 
boundary conditions, along the compactified Euclidean time direction. Then ao's become the 
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coordinates of D-instantons on the dual circle. The integration over ao's becomes now the 
integration over the positions of D-instantons. The partition function should involve such 
an integration over the collective coordinates and since they are collective coordinates the 
integration appears in front of the exponential of the effective action, not in the action itself. 
Viewed in terms of DO-branes and open strings, we have a gas of DO-branes with open strings 
between them. The individual strings might have a winding number q (more precisely 2q + l 
in the case of superstrings), describing the winding around the finite-temperature space-time 
cylinder. The energy of such states are oc j3q/2na'. However, the g's satisfy J2q — as a 
result of the integration over ao- Physically this constraint is most easily understood by 
going to the closed string channel where we have closed string boundary state on the dual 
circle with radius (3 = ATt 2 a'/f3 localized at the point (u$,q). Passing to the momentum 
representation, we write 



same energy as the above mentioned open string states. In this representation J2 Q — 
simply expresses momentum conservation in the thermal direction. 

The effective static potential between two DO-branes emerges because supersymmetry is 
broken by finite temperature. This effect of breaking supersymmetry is somewhat analogous 
to the velocity effects at zero temperature where the matrix theory and superstring compu- 
tations agree to the leading order of the velocity expansion!. We have thus shown that the 
leading term in a low temperature expansion is correctly reproduced by the matrix theory. 
The discrepancy between the matrix theory and superstring computations, which we have 
observed in the limit of large distances LT ^> 1, does not contradict this statement since 
temperature the limits of large distances and small temperatures are not interchangeable. 

An interesting feature of the effective static potential between DO-branes is that it is 
logarithmic and attractive at short distances. In the matrix theory, the singularity of the 
computed 1-loop potential occurs when the distance between the DO-branes vanishes and 
the SU(N) symmetry which is broken by finite distances is restored. The integration over 
the off-diagonal components can no longer be treated in the 1-loop approximation! In the 
superstring theory, the singularity is exactly the same as in the matrix theory since it is 
determined only by the massless bosonic modes (the NS sector in the superstring theory). 
Its origin is not due to the presence of massless photon states in the spectrum. Putting 
v\ = v-i in the above D-instanton picture on the dual temporal circle, we see that the mass 
of the lowest states, associated with the winding numbers 2q + 1 = ±1 is f3/2ira' . The 
divergence at L = emerges, in this picture, after summing over all the open string states 



oo 




(45) 



Here the temporal momentum is quantized as po = 2irq//3 = qj3/2Tca', which lead to the 
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since no single state has such a divergence. It shows up only at finite temperature where the 
winding number q exists. 

It is important to notice that the computed partition functions take into account only 
thermal fluctuations of superstring stretched between DO-branes but not the fluctuations of 
DO-branes themselves. To calculate the thermal partition function of DO-branes, a further 
path integration over their periodic trajectories a(r) is to be performed as in (^). One might 
think that classical statistics is applicable to this problem since the DO-branes are very heavy 
as g\ M — >• so that one could restrict himself by the static approximation. This is however 
not the case due to the singularity of the effective static potential at small distances. The 
integral over the DO-brane positions efs is divergent when the two positions coincide. 

However, this singularity is only in the classical partition function. The path integral over 
the periodic trajectories a(r) that we actually have to do can not diverge since the 2-body 
quantum mechanical problem has a well-defined spectrum. The path integral can then be 
evaluated as J2n ex P( — PE n ) where E n are in the spectrum of the operator H = P 2 / M + V e s. 
There certainly should not be the bound state energy eigenvalue at negative infinity for this 
quantum mechanical problem which implies the convergence of the path integral. These 
issues which are related to thermodynamics of DO-branes will be considered in a separate 
publication. 

Let us finally discuss when the 1-loop appoximation that we have done is applicable. The 
loop expansion in Yang-Mills theory computation is valid only in the limit where 



is small. This is due to the fact that the distance L plays the role of a Higgs mass which 
cuts off the infrared divergences of the loop expansion in the 0+1 -dimensional gauge theory. 
Thus, the perturbative Yang-Mills theory computation is good when 



This can be satisfied if either the string coupling is small or if the DO-brane separation is large 
compared with the string length scale. In the latter case, the truncation of the spectrum to 
the lightest modes is still valid when 



Note that the first inequality which is independent of both the string scale and the string 
coupling is the one already discussed in the previous section. In this case, the temperature 
must be less than the inverse distance between DO-branes. In the case where the string 





(46) 



k B T < - < 



1 



(47) 
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coupling g s is small, the criterion for validity of truncation of the spectrum becomes 



k B T < 



(48) 



that is the usual one. 
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